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Abstract 

The probability distribution of random walks on linear structures generated 
by random walks in d-dimensional space, P d (r, t), is analytically studied for 
the case £ = r/i 1 / 4 <C 1. It is shown to obey the scaling form P d (r,t) = 
/?(r)t _1 / 2 ^ _2 /rf(^), where p(r) ~ r 2 ~ d i s the density of the chain. Expanding 
fd(0 m powers of £, we find that there exists an infinite hierarchy of critical 
dimensions, d c = 2,6, 10,..., each one characterized by a logarithmic correc- 
tion in f d {£). Namely, for d = 2, / 2 (£) ~ a 2 £ 2 ln£ + 6 2 £ 2 ; for 3 < d < 5, 
fd(0 - ad 2 + b d £ d ; for d = 6, / 6 (£) ~ a 6 £ 2 + & 6 £ 6 ln£; for 7 < d < 9, 
fd(0 -a d e + b d f + c d £ d ; W= 10, f 10 (£) ~ ai £ 2 + &io£ 6 + cio£ 10 lnf, etc. 
In particular, for d = 2, this implies that the temporal dependence of the prob- 
ability density of being close to the origin Q2(r, t) = P2(r, t) / p(r) ~ i -1 / 2 Int. 

Pacs: 5.40.+j, 05.60. +w, 66.30.-h 
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I. INTRODUCTION 



Random fractals represent useful models for a variety of disordered systems found in 
Nature. In addition to their structural properties, fractals have attracted much attention in 
recent years because of their interesting transport properties £[]-|J . 

Of particular interest is the question of how the probability density of random walks, 
P d (r,t), is changed on fractal structures with respect to its Gaussian form valid on regular 
d— dimensional systems, P d (r,t) ~ t~ d l 2 exp(— const x r] 2 ), where 77 = r/t 1 ^ 2 . The form of 
Pd(r,t) on fractals has been extensively studied in the asymptotic limit £ = r/t l ^ dw 3> 1 
PJ5|-[T2|, where d w is the anomalous diffusion exponent characterizing the time behavior of 
the random walks, (r 2 (t)) ~ £ 2 M". As a result of these investigations, it is now generally 
accepted that P d (r,t) displays a stretched Gaussian form 

P d (r,t) ~ p(r)t- da/2 exp(-const x £"), £ > 1, (1) 

where p(r) ~ r dj-d - g density of the fractal structure, df is the fractal dimension, 
d s = 2df/d w is the spectral dimension u = d w /(d w — 1), and is normalized according 
to / dr r d ~ 1 P d (r,t) = 1. However, much less is known about the behavior of P d (r,t) in the 
opposite limit when £ approaches zero. 

In this paper we concentrate on diffusion in linear random fractal structures generated 
by random walks (random-walk chains, RWC) in d- dimensional systems, where Pd(r,t) can 
be obtained exactly. Recently, using numerical simulations, it has been suggested that for 
such linear fractals [|l"3|j , 

Qd(r, t)/Q d (0, t) ~ (1 - const x ^~ 2 ), £ - 0, (2) 

for all dimensions d, where Qd{ r ,t) — p( r )~ 1 -Pd( r ) t) is normalized on the fractal chain, 
i.e., J drr^- l Q d (r,t) = 1, with d f = 2, d w = 4 for RWC, £ = r/i 1 / 4 and Q d (0,t) is the 
probability density to return to the origin. 

In the following, we derive an exact expansion for P d {r, t) in the limit of £ — > 0. Surpris- 
ingly, P d (r, t) displays an extremely rich behavior as a function of both £ and dimensionality 
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d. We show, among other results, that Eq. (0) can only be valid for 3 < d < 5, and 



P d (r,t) ~ p(r)t 



/2 (1 -const x £ 4 ), for d > 7. 



(3) 



Moreover, we find that the small-£ expansion of Pd{r, t) is characterized by a hierarchy of cri- 
tical dimensions, d c = 2, 6, 10, 14, . . ., where logarithmic corrections of the form £ dc ~ 2 log(l/£) 
occur. In particular, for d = 2 we obtain P 2 (r,t) ~ 2p(r)t -1 / 2 ln(t 1 / 4 /r). 



We consider linear structures generated by random walks in <i-dimensional systems. Such 
structures are fractals with fractal dimension df = 2, independently of d. To study diffusion 
of particles along such linear chains, we assume that the diffusing particles (random walk- 
ers) can move only along the structure (path) which has been created sequentially by the 
generating walks. Thus, although the structure can intersect itself in space, the walkers see 
just a linear path. We denote such paths as random- walk chains (RWC). 

Along the linear path, the probability density of random walkers, at chemical distance £ 
along the RWC from their starting point after time t, p(£, t), subject to the initial condition 
p(£,0) = S(£), approaches the well-known Gaussian distribution 



normalized according to / °° d£p(£,t) = 1. Thus, diffusion along the chain (i.e., £-space) is 
normal and (£ 2 ) = t. On the contrary, in Euclidean r-space diffusion is anomalous with 
d w = 2d f = 4 (see, e.g., @). 

To obtain the behavior of the probability density in r-space, averaged over all RWC 
configurations, Pd(r,t), we note that it is related to p(£,t) by 



and is normalized according to / d d r Pd(r,t) = 1. Another possibility is a normalization on 
the RWC fractal, i.e., dr r d f~ l Qd(r,t) = 1. Both distributions are simply related to each 



II. RANDOM WALKS ON RANDOM- WALK CHAINS 




(4) 




(5) 



other by P d (r,t) = p(r)Q d (r,t). 



3 



In Eq. (H), $rf(r, £) represents the probability for a site r to belong to a RWC at distance 
I from the origin along the chain. The chemical distance i plays the role of the time variable 
in Eq. (II), and one can immediately write 



A n 



1 \d/2 



2nt 



6XP [-21. 



(6) 



where is a normalization factor such that / d d r Q d (r,£) = 1. Therefore, by inserting 
and (|) in (g) we infer @,||] 

1 \ d / 2 2A d /•«> / r 



2W (2^o Mrd/2 ^(-Q^(4t 



(7) 



Now, the elementary transformation x 



P d (r,t) = 2A d (2ir) 



brings (|7[) to the form 



-(d+l)/2 -d 



(8) 



where the scaling function /<i(£) is defined by 

m)=? dxx^exp 
Jo 



(9) 



for the scaling variable £ = r/t 1 / 4 . If the RWC normalization is chosen, the distribution 
Q d (r,t) = p-\r)P d {r,t) r 1 /^), where / d (£) = r 2 /d(0- 

To deal now with the evaluation of f d {C) when £ — > 0, it is convenient to rewrite the 
integrand exponent as 



exp 



-- Uv + I 

2 V x 



exp 



exp 



2x V a; 



and expand the second exponential factor in Taylor series. The remaining integrals can be 
solved exactly (see, e.g. fI5fl), and one arrives at the following expression for (0) 



oo -| ✓ i x n n 

/-co = e e ' 



n=0 ' t - V Z/ fc=0 



W Ud/2-l 



E(-i) fc U]e d/2 - 1+ ^ W2 -i + . +fc )(e 2 ), 



(10) 



where if v is the modified Bessel function of order v. 

Let us consider Eq. flT0p in some particular cases of interest. The results for spatial 
dimensions d < 7 are summarized in Table |. All the coefficients were calculated numerically 



by computing the double sums explicitly. In some cases they are available in analytic form, 
but we include their numerical values to make the table uniform. 

However, besides the coefficients, the main properties of these expansions can be obtained 
readily as follows. The key parameter is s — |(| — 1). The corresponding values of d = 
2(2s + 1) for integer s should be referred to as critical dimensions, d c = 2, 6, 10, . . .. Each 
order in the expansion has its own critical dimension. The leading term has d c = 2, the 
first correction term has d c = 6, the second correction term has d c = 10, etc. This has to 
do with the functional form of fd(£) in the corresponding order which for d < d c depends 
on d, at d — d c it has a logarithmic correction and for d > d c becomes independent of d. In 
particular, the leading term of fd(£,) behaves as £ for d = 1 and £ 2 ln(l/£) for d — d c — 2 
and as £ 2 for all d > d c = 2, the first correction term behaves as £ d for 2 < d < 6, and for 
d = d c = 6 as £ 6 ln(l/£) and as £ 6 , for all d > d c = 6, and so on. 

Mathematically, this behavior can be explained by the intrinsic properties of the Bessel 
function K s (£ 2 ). By its definition, K s (^ 2 ) = |cosec(7rs)[/ s (£ 2 ) — J_ s (£ 2 )] for non-integer s 
and, in turn, I s (^ 2 ) = £ 2s ET=o b k(s)^ k and £ 2s /- s (£ 2 ) = ET=obk(- s )^ 4k ■ A s one can see 
this expansion has s-independent powers of £ that form the invariant part of fd(0- The 
first terms of £ 2s K s (l; 2 ) expansion are 

e s K s (e) « b (s)t + h(s)t +i + ■■■- bo{- 8 ) - hi-s)^ 

Thus, for < s < 1 (i.e., 2 < d < 6) f d (g) has the form / d (£) w ^ 2 [a (-d) - a {d)£ d - 2 ). We 
see that the first term of this expansion is the invariant part (up to numerical coefficients) 
of which remains unchanged when varying d. The same argument shows that for 

1 < s < 2 (6 < d < 10) / d (0 takes the form f d (£) w £ 2 {a (-d) - a x (-d)£ 4 + a (d)£ d - 2 } and 
now two first terms of this expansion are the invariant part of f d (£). A special case in our 
problem is d — 1, i.e., s — — \ < 0. Then the leading term of f d (£) is £ 2 £ -4 l s l = £, which is 
easily seen by noting that _fT_ s (£ 2 ) = K s (^ 2 ). 

For integer values of s, the small-£ expansion of £ 2s K s (!; 2 ) has a logarithmic term of form 
e 4s log(l/0=^- 2 log(l/0 0. 
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In general, there are [s] + 1 ([s] is the integer part of s) terms in the invariant part of 

U0- 



III. CONCLUSIONS 

We have studied analytically the small-£ expansion of the mean probability density, 
P d (r,t), of random walks on random-walk chains in ci-dimensional space. We have shown 
that the leading terms of the expansion of, P d (r,t), behaves, in the limit £ = r/t 1 / 4 — > 0, as 

P d (r, t) oc p(r)r 1/2 (l - a^ 2 ), when 3 < d < 5, 

and as 

P d (r, t) oc p(r)r J / 2 (l - q£ 4 ), when d > 7, 

where p(r) ~ r d /- rf anc [ ^ — 2. This implies that the probability density Q d (r,t) = 
P d (r,t)/p(r) on the fractal chain behaves for d > 7 as Q d (r,t) ~ t _1 / 2 (l — c d r 4 /t), consistent 
with the behavior of diffusion in £-space, i.e., p(£,t) ~ t _1 / 2 (l — £ 2 /2t), for I 2 <^ t, and the 
fact that Q d {r — ~ = when d — > oo. We see that this already occurs when 
d > 7. 

We have shown that logarithmic corrections occur at critical dimensions d = d c = 4n + 2, 
with n — 0, 1, 2, . . ., i.e., d c = 2, 6, 10, . . ., for the terms £ dc ~ 2 log(l/£). In particular for d = 2, 
Q 2 (r : t) ~ i -1 / 2 ln(t 1/,4 /r), for r <^ t 1 / 4 , and the probability density for the random walker 
to be close to the origin, Q 2 (r,t) behaves as t _1 / 2 lnt. This logarithmic correction is due to 
the fact that in two dimensions the RWC returns to its starting point with probability one. 
In one dimension, Q\(r, t) ~ i -1 / 4 /r, and in d = 3, Q 3 (0,t) ~ t -1 / 2 . One can say that d — 2 
plays the role of a marginal dimensionality for the probability density of being at the origin 
of random walks on RWC, while for larger r each order in the expansion has its own critical 
dimension. 
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TABLES 

TABLE I. The leading and correction terms for the series expansion of fd(Oi with £ 
when £-> 0, Eq. (Mj, as a function of dimension d. 



= r/t 1 ^, 



d Leading term First correction Second correction 

1 2.1558 i -2.5066 £ 2 0(£ 5 ) 

2 2£ 2 ln(l/0 0.1738£ 2 0(£ 6 ) 

3 2.5066 £ 2 -0.0609 £ 3 0(£ 6 ) 

4 2 £ 2 -1.2533 £ 4 0(£ 6 ) 

5 2.5066 £ 2 -1.4372 £ 5 0(£ 6 ) 

6 4£ 2 -£ 6 ln(l/£) -0.3369£ 6 

7 7.5199 £ 2 -1.2533 £ 6 0.8244 £ 7 
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